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by Fred Teren and Gary L. Cole 
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SUMMARY 

This report describes a simplified analytical method for calculating partial deriva- 
tives relating midcourse correction velocity requirements to injection errors resulting 
from guidance system errors for lunar and planetary missions. 

The analytical method uses two-body equations of motion to describe the reference 
transfer orbit. These equations are linearized to obtain partial derivatives relating 
changes in terminal state conditions to changes at injection and midcourse. The partial 
derivatives are then used to calculate the required correction velocity. 

The analytical equations are easily programmed on a digital computer and require 
very little computer time. The only input variables required are the position and ve- 
locity components at injection from a reference trajectory. 

The analytical partial derivatives are compared to exact partial derivatives obtained 
from a detailed N Body computer simulation. A comparison of results is presented for 
two lunar and five planetary missions. The error in the approximate results is found to 
be less than 14 percent for all cases simulated. 


INTRODUCTION 

The predicted guidance system performance for a space mission is usually mea- 
sured in terms of a single number, called the figure of merit (FOM). For lunar and 
planetary probes, the FOM is the expected value of the midcourse velocity correction, 
which is used to correct miss plus time of flight or miss only errors at the target. 

The FOM may be calculated as the product of two matrices: the injection covariance 
matrix (ICM) and the FOM covariance matrix. The ICM gives a complete statistical 



description of the injection errors and may be obtained by integration of the adjoint of the 
equations of motion from launch to termination of powered flight for each guidance system 
error source. The FOM covariance matrix consists of a set of partial derivatives re- 
lating FOM to the injection errors. The exact calculation of these partial derivatives is 
a tedious one requiring an N Body computer program and much computer time. The pur- 
pose of this report is to derive an approximate analytical method which allows these par- 
tial derivatives and, hence the FOM covariance matrix, to be easily calculated. 

In preliminary design studies, the mission planner may wish to investigate the FOM 
requirements for various combinations of trajectories, launch vehicles and/or guidance 
hardware error models. In such cases, the detailed mission requirements and launch 
opportunities usually have not been established; and targeted reference trajectories are 
not available. Also, the mission planner may not have a detailed N Body code at his dis- 
posal, and therefore could not calculate exact (N Body) partial derivatives. It is for this 
type of problem that the authors decided to investigate possible analytical methods for 
generating the partial derivatives required to compute the FOM. The analytical methods 
would be easy to implement on a computer and require little computer time. 

As a result of this study, equations were developed for computing analytical partial 
derivatives of FOM with respect to injection state variables for planetary and lunar 
probes. The analytical equations for computing these derivatives are the same for the 
two mission types, but the approach used is somewhat different. In both cases, changes 
in arrival conditions at the target are related to injection errors and to midcourse cor- 
rections by using a method described by Danby in references 1 and 2. Basically, 

Danby's method relates errors between any two points on a two -body orbit by using a 
first-order analysis. 

A digital computer program was written to test the analytical equations and to com- 
pare the analytical results to the detailed N Body results. A reference powered flight 
trajectory is required as input to the program. In addition, a typical guidance hardware 
error model and error values were assumed. These errors provide the forcing function 
for the linearized equations of motion, which are integrated in the program to obtain the 
injection covariance matrix. The reference transfer orbit is obtained by using the injec- 
tion conditions from the reference powered flight trajectory in the two-body analytical 
equations. 

Analytical and N Body results are presented for lunar and planetary missions, in- 
cluding Venus, Mars, Mercury, and Jupiter. The results presented include the analyt- 
ical and N Body FOM for each error source, as well as the total FOM. The analytical 
and N Body FOM covariance matrices are also presented. A sample problem is pre- 
sented wherein the trajectory resulting in the largest FOM for a Mars mission is deter- 
mined by using the analytical technique. This problem serves to illustrate the use of the 
analytical technique. 
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ANALYSIS 


General Equations for FOM 

For planetary and lunar probe missions, the function of the launch vehicle guidance 
system is to steer the vehicle during its powered phases and to command cutoff, such 
that the resulting free flight trajectory will arrive at the desired target, usually at some 
prespecified time. The guidance system performance for these missions can be specified 
in terms of a single number, called the figure of merit (FOM). The FOM is the expected 
value of the midcourse correction velocity, required to correct miss plus time of flight 
errors at the target body or to correct miss only, if time of flight errors are not impor- 
tant. 

For a definite set of injection errors, the required midcourse correction can be de- 
termined by simulating the dispersed trajectory to the target, and iterating on the mid- 
course correction to zero the error in arrival conditions. Since injection position and 
velocity errors due to guidance are usually small, a first-order analysis has generally 
been used to propagate state vector deviations from injection and midcourse to the target. 
The use of first-order equations also allows the effects of various guidance system error 
sources to be statistically combined to obtain the FOM. 

If a first-order analysis is used, the midcourse correction for a set of injection er- 
rors can be expressed as 


Av c = G5S 

where G is a 3 x 6 matrix of partial derivatives 

a ( Av c), i = 1, 2, 3 

G ik = 

as k k = 1 , . . . , 6 

and 


5S = 



Injection is defined here as the time of termination of powered flight. 
The correction velocity squared can be expressed as: 

(Av c ) 2 = AvJ Av c = (G6 s) T (G6S) 


( 1 ) 


( 2 ) 


3 



where T denotes the transpose of a matrix. Equation (2) may be rewritten as: 


(Av c ) 2 = 6S T G T GSS (3) 

T 

The 6x6 matrix G G is defined as the FOM covariance matrix 

A = G T G 
T 

By replacing G G with A, equation (3) may be rewritten in summation form as: 

6 6 

K) 2 = Z Z v s i 6s i (4 > 

i=l i=i 

The expected value operator E can be applied to equation (4) to give: 

(FOM) 2 = E [(4V C ) 2 ] = £ 2 A iflj (5) 

i=l j=l 

where 

cT ±j = E^SjSSj) 

The same A matrix applies for either a definite or a statistical set of injection errors. 


Calculation of FOM Covariance Matrix - Planetary Mission 

For a planetary mission, the FOM is the midcourse velocity required to correct 
miss plus time of flight errors at the target planet. The correction is usually performed 
several days after injection, after enough tracking data has been obtained to define the 
transfer orbit. 

For the analytical calculations, a series of two-body transfer conics is assumed, 
with the transition from the earth's field to heliocentric space occurring at the earth's 
sphere of influence. The midcourse correction is assumed to take place at the earth's 
sphere of influence. This gives a correction time of from 1 to 4 days (depending on the 
injection energy) which is consistent with the real case. The following assumptions are 
also made: 
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(1) The error in time at earth escape is negligible compared to the total trip time. 

(2) The position error at earth escape is negligible compared to heliocentric dimen- 
sions. 

With these assumptions, the effect of injection position and velocity errors is simply to 
change the magnitude and direction of the hyperbolic velocity vector. Therefore, the 
midcourse correction velocity is equal to the error in the hyperbolic velocity vector. 

The first step in the analysis is to compute the reference trajectory parameters at 
injection. All reference orbit parameters in this report are calculated from the two- 
body equations tabulated in reference 3. Some of the orbit parameters defined in equa- 
tion (6) are illustrated in figure 1. 



Figure 1. - Trajectory geometry for Earth-centered phase of planetary mission. 
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Angular momentum: 


*\ 


Similatus rectum: 


Energy: 


Eccentricity: 


h = r x v 



E = XJLY. - i£ 

2 r 


Perigee radius: 


Injection true anomaly: 


Asymptote true anomaly: 


Time from perigee: 



> ( 6 ) 


t-t p = 'Xl — 




!og 


' Ve + 1 + Ve - 1 tan - 

2 1 e sin r ? 


V^Ttan^/ 1 + 6 cos t] 
2 ; 


All symbols are defined in appendix A. 

The error in the hyperbolic velocity vector is calculated by Danby's method as de- 
scribed in references 1 and 2. Danby’s method is briefly described in this report in 
appendix B. 
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For the planetary case, interest lies in calculating the 6v at earth escape. 


Av c = 6v e = C eI 6r I + D eI 6v I 


where C g j = C(t e ,tj) and D gI = D(t g ,tj) are calculated from equation (B2) in appendix B. 

The form of equation (1) has been satisfied by equation (7) by combining C and D 
into the correction sensitivity matrix G: 

G = (C | D) 

The A matrix is then given by 

A = G T G 


By examining the equations in appendix B, it will be noted that the L and M ma- 
trices go to infinity as time approaches infinity. However, the P and Q matrices re- 
main finite. Physically, this means that the position errors increase without bound for 
any injection errors. However, in keeping with the assumptions made earlier, the posi- 
tion errors are assumed negligible with respect to heliocentric dimensions and assumed 
to be zero for this case. This will be a good assumption if the hyperbolic velocity is 
large so that earth escape occurs quickly and position errors do not have time to build 
up appreciably. 

It is necessary to compute the unit radial coordinate system at perigee since Danby’s 
equations are referenced to perigee. This is done as follows: 

r = 1 - ** - JL 
P lie er 



A A A 

W r P 

where r, v, and h are determined at injection. All state vector deviations must be 
expressed in this coordinate system. 

The complete set of equations required to calculate the A matrix for the planetary 
case is presented in appendix C. 
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Calculation of FOM Covariance Matrix - Lunar Mission 

The FOM for the lunar case is the expected value of the midcourse velocity required 
to correct either miss plus time of flight to the moon or miss only. This correction is 
usually made about 15 to 20 hours after injection, but can be computed for any desired 
correction time by using either the exact or analytical methods. 

The Av c is calculated by using Danby's matrices - the same technique that is used 
in the planetary case and described in appendix B. The miss at the moon is calculated 
at the incoming asymptote at the lunar sphere of influence, relative to the reference 
geocentric conic trajectory. In this case, the following definitions are made: 

t = time of midcourse correction 
t^ = time of lunar arrival 
The uncorrected miss at the moon is given by: 

6T m = Aml^I + B mI 6 ^I + < ? A " V 6t A 

where A mI = A(t m ,tj) and B mI = B(t m ,tj) are calculated from equation (B2) in appen- 
dix B. 


v^ = arrival velocity at the moon 


v m = moon's orbital velocity about the earth 


6t^ = error in arrival time 


The change in arrival position due to midcourse correction is: 


6r = B 6v 
m me c 


To correct miss plus time, the time error 6t^ must be zero. Therefore, 


A? c, MT = - B mc ( A m^I + B ml^l) = G6S s U 


(8) 


8 



where 


G =-l B mc A mI 


B" 1 B T 
me ml 


Therefore, 


A MT ' G g 


To correct miss only, the arrival time t A may be optimized to give minimum Av f 

'c,M 


A V M ■ U - B^ c ( 7 A - *m>«A = » + V5t, 


where 


V = -B” ^ ( v » - v ) 
me' A nr 


The following is obtained by squaring equation (9): 


For optimum t A , 


Therefore, 


( Av c, m) 2 = uTu + 2U T V6t A + V T V6t A 


8 ( At c,m ) 2 _ 

at a 


= 2(u r v + v^vst^ = o 


5t , 


ITV 

V T V 


and 


ITi — \2 




V T V 


The A matrix can be constructed for the miss only case by using equation (8): 



1 1 I It III 


III III I I 


i an ii n ii i lain iii i in hi i 


i i ii i 



<^ T 4 s 

v 1 ^ / 


where G refers to the miss plus time matrix; therefore, from equation (3), 


A M " A MT 


G^V V T G 


V"V 


for the miss only case. 

As in the planetary case, the elements of the A and B matrices are calculated 
from the equations in appendix B. The reference trajectory parameters are the same as 
for the planetary case (eq. (6)), except that true anomaly at lunar arrival is used in place 
of true anomaly at the outgoing asymptote. 

True anomaly at lunar arrival: 


,-l 




The time from perigee equation is: 


t 




1 - e 
1 + e 



e sin r] 

1 + e cos r\ 


Some of the trajectory parameters are illustrated in figure 2. For the lunar mission, 

v» and v_ are calculated as follows: 

A m 

Unit vector to the moon: 


r A - cos Vp + sin Vp 


Velocity at lunar arrival: 


v. = - h X r 
A p P 
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To calculate the lunar velocity, the following are assumed: 

(1) Lunar orbit is circular about the earth. 

(2) if/ is the inclination of the moon’s orbit to the earth’s equator. 
As shown in appendix D, these assumptions result in 


V rv, = 

m 


■,/Z 1 f cost a xi) ± -i/l r; x( ; x S) ]l 

y r m co® „ | cos „ - y cos^<p m m J j 


where w is the latitude of the moon at arrival: 


sin <p = r m • z 


The choice of sign in v m depends on whether the moon is ascending or descending: 
+ sign indicates descending moon 
- sign indicates ascending moon 
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The complete set of equations required to calculate the A matrix is presented in appen- 
dix C. 


Calculation of Injection Errors 

In order to test the analytical equations for FOM, it was necessary to generate an 
injection covariance matrix. This was done by choosing a typical guidance hardware 
error model and integrating the linearized equations of motion. 

Guidance hardware error model. - A typical platform inertial guidance system 
hardware error model was chosen and is described in table I. Among the errors con- 
sidered were accelerometer bias and scale factor, gyro bias and g-sensitive drift and 
platform misalinement. The errors, when applied to the nominal thrust acceleration 
profile, provide the forcing function 6a for the linearized equations of motion. The 
platform is alined along the u, v, w inertial coordinate system which is established at 
liftoff. The u axis points in the pad azimuth direction. The w axis lies along the local 
plumbline and v completes the right-handed orthogonal set. 

Linearized equations of motion. - The equations of motion for the reference trajec- 
tory in vector form are: 


r = 



(II) 


— „ 3 __ 

where r is the total acceleration, -11 r/r is the gravity acceleration, and a is the 

thrust acceleration. If the error in sensed acceleration 6a is small, then the error in 
total acceleration (i. e. , difference between measured and actual acceleration) may be ob- 
tained by linearizing equation (11). The result is: 



6r + — (r • 5r)r + 6a 


(12) 


The difference between the measured and actual state vectors, 8S, is obtained by inte- 
grating equation (12). 

Injection error s. - The injection state vector error is defined by 

rf-Acf Vl = - 5 ^<WS ref ) (13) 

where the subscripts (act, ref, mes) refer to the (actual, reference, measured by the 
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TABLE I. - DESCRIPTION OF HARDWARE ERROR MODEL 


Error source 

Error 

symbol 

Description 

Error 

value 

Error source 

Error 

symbol 

Description 

Error 

value 

Initial platform misalinement, 

e i 

u Axis 

0.054 

Absolute accelerometer 

e 25 

u Accelerometer 

. 0 

mrad 

*2 

v Axis 

.054 

scale factor, percent 

e 26 

v Accelerometer 

0 


*3 

w Axis 

.090 


e 27 

w Accelerometer 

0 

Accelerometer input axis 

e 4 

v Relative to u 

0.050 

Accelerometer scale factor 

e 28 

u Accelerometer 

0.0094 

misalinement, mrad 

e 5 

w Relative to u 

.050 

(quadratic term), mg/g^ 

e 29 

v Accelerometer 

.0094 


e 6 

w Relative to v 

.055 


e 30 

w Accelerometer 

.0094 

Accelerometer preflight 

e 7 

u Accelerometer 

0.042 

Accelerometer scale factor 

e 31 

u Accelerometer 

0 

bias error, mg 

e 8 

v Accelerometer 

.042 

(cubic term), mg/g 3 

e 32 

v Accelerometer 

0 


e 9 

w Accelerometer 

.042 


e 33 

w Accelerometer 

0 

Accelerometer scale factor 

e 10 

u Accelerometer 

0.051 

Accelerometer cross-coupling 

e 34 

u Accelerometer 

0.013 

(linear), percent 

e ll 

v Accelerometer 

.051 

(input -pendulous axes), mg/g 2 

e 35 

v Accelerometer 

.013 


e 12 

w Accelerometer 

.051 


e 36 

w Accelerometer 

.013 

Gyro fixed torque drift, 

e 13 

u Gyro 

0.084 

Accelerometer cross-coupling 

e 37 

u Accelerometer 

0.012 

deg/hr 

e 14 

v Gyro 

.093 

(input -output axes), mg/g 2 

e 38 

v Accelerometer 

.012 


e 15 

w Gyro 

.094 


e 39 

w Accelerometer 

.012 

Spin axis mass unbalance 

e 16 

u Gyro 

0. 173 

Accelerometer inflight bias 

e 40 

u Accelerometer 

0.024 

drifts, deg/(hr)(g) 

e 17 

v Gyro 

. 190 

error, mg 

e 41 

v Accelerometer 

.026 


e 18 

w Gyro 

. 177 


e 42 

w Accelerometer 

.029 

Input axis mass unbalance 

e 19 

u Gyro 

0. 106 

Accelerometer bias error 

e 43 

u Accelerometer 

0.009 

drifts, deg/(hr)(g) 

e 20 

v Gyro 

.101 

output axis loading (quadratic), 

e 44 

v Accelerometer 

.009 


e 21 

w Gyro 

. 114 

mg/g 2 

e 45 

w Accelerometer 

.009 

Anisoelastic drifts, 

e 22 

u Gyro 

mm 

' Accelerometer cross -coupling 


u Accelerometer 

1 SB 

deg/(hr)(g 2 ) 

e 23 

v Gyro 


(pendulous -output axes), 


v Accelerometer 

.008 


e 24 

w Gyro 


mg/g 2 


w Accelerometer 

.008 









guidance system) trajectories. If a guidance cutoff is used and perfect guidance software 
is assumed, then the midcourse correction requirement results solely from hardware 
errors, and 


Av =GAS= -G6S 

Therefore, the midcourse correction may be obtained by applying G to the 6S obtained 
by integration of equation (12). A guidance cutoff will normally be used for the lunar and 
planetary missions in order to obtain the desired injection energy. 


Calculation of Injection Covariance Matrix 

The 6x6 matrix a is called the injection covariance matrix (ICM). The injection 
errors result from errors in the guidance system accelerometers and gyros. Since the 
expected hardware errors are of a statistical nature, the elements of the ICM are defined 
as the expected values of the products of the state variable deviations; that is: 

a ti = ECSSjSS.) i = 1, . . • , 6 

1J 13 j = 1, • • • , 6 

The state variables are the reference trajectory position and velocity components at 
injection and the 6S^ are the corresponding deviations in these state variables. 

For each guidance system error source, the injection state deviation can be ex- 
pressed as: 


5S i " P ik e k 


i = 1, 


6 


1U 

where e k is the magnitude of the k in error and P ik is the matrix of partial deriva- 
tives 9S i; /3e k . The elements of the injection covariance matrix (for N guidance system 
error sources) are given by: 


N N 

a - EfaSjSSj) = £ 2 
k=l 1=1 

If the various error sources are assumed to be independent, then the elements of the 
injection covariance matrix become: 
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(14) 


N 

“ii = E < p ik p )k) E (4) 


k=l 


2 

The total (FOM) is obtained by combining equations (5) and (14). 

(FOM) 2 = E(eI A P ik p\(e 2 ) (15) 

k=l\i=lj=l / 

For any given error model, the quantities in parentheses in equations (14) and (15) rep- 
resent the contribution of each error source to the injection covariance matrix and to 

2 2 
(FOM) , respectively. Once these coefficients are calculated, the ICM and (FOM) can 

be obtained for any error budget by performing the indicated summations in equations 
(14) and (15). The contribution of the individual error sources can also be easily calcu- 
lated. 

In the computer program, the linearized equations of motion are integrated to obtain 
6S for each error source. A FOM is then computed for each error source. The total 
FOM is obtained by calculating the RSS of the individual FOM's. 


RESULTS AND DISCUSSION 

In order to determine the accuracy of the approximate analytical equations pre- 
sented, both analytical and exact (N Body) error analysis results have been obtained and 
are presented here for comparison. For the N Body results, targeted powered flight and 
free flight trajectories were obtained by using the procedures discussed earlier. The 
lunar trajectories were targeted for a vertical impact with a prespecified time of flight. 
The planetary cases were also targeted to a prespecified flight time, but a close plane- 
tary approach was used as the targeting criteria. After the targeted reference trajectory 
was obtained, the linearized equations of motion were integrated from powered flight in- 
jection to the target body and from the midcourse correction point to the target. The 
partial derivatives obtained in this manner were then combined to give the exact FOM 
covariance matrix, as discussed earlier. 

The analytical FOM covariance matrix is obtained by using the injection conditions 
from the targeted powered flight reference trajectory in the analytical equations pre- 
sented earlier. It is important to note that injection conditions could have been obtained 
from many sources, such as references 4 and 5. As stated earlier, it is not necessary 
to have targeted reference trajectories in order to obtain analytical error analysis re- 
sults. The targeted injection conditions are used here to insure an exact comparison 
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between the analytical and N Body results. 

The position, velocity, and acceleration profiles from the reference trajectory are 
required as input to the error analysis program. Several different configurations of the 
Atlas-Centaur launch vehicle were used to simulate the powered flight reference trajec- 
tories, as shown in table n. These configurations do not necessarily represent existing 
or planned launch vehicles or stages; however, the configurations used are realistic 
ones, and the results obtained can be considered representative. 


TABLE B. - GENERAL DESCRIPTION OF TRAJECTORIES SIMULATED 


Case 

Mission 

Ascent 

mode 

Launch 

date 

Launch 

azimuth, 

deg 

1 

Mars 

Direct 

3/22/69 

111 

2 

Mars 

Direct 

12/13/66 

120 

3 

Jupiter 

Direct 

2/26/72 

95 

4 

Venus 

Two -burn 

8/13/70 

90 

5 

Mercury 

Two -burn 

11/12/68 

90 

6 

Lunar 

Direct 

6/3/66 

93 

7 

Lunar 

Two -burn 

1/1/69 

114 


a High-energy cryogenic upper stage. 


Energy, 

km 2 /sec 2 

True 

anomaly, 

deg 

Parking 

orbit 

coast 

time, 

sec 

Launch vehicle 

14.45 

-7.54 

— 

At las -Centaur - 
upper stage a 

15.07 

16.67 

— 

Atlas -Centaur - 
upper stage a 

114 

-2.41 

— 

Atlas-Centaur - 
Burner n 

15.05 


445 

Atlas -Centaur 

45.36 


1232 

Atlas -Centaur 

-1.20 

.03 

— 

Atlas-Centaur - 
upper stage 3. 

-1.22 


1201 

Atlas-Centaur 


For the two-burn lunar, Venus and Mercury missions, an Atlas Centaur vehicle 
was simulated. The powered flight trajectory started with a boost into a 100-nautical 
mile parking orbit. The Atlas and Centaur were used to attain orbital altitude and 
velocity. The vehicle then coasted in orbit until the vehicle had ’’caught up” with the 
desired transfer conic. At this time, the Centaur stage was reignited and burned until 
the required final injection conditions had been achieved. Steering of the launch vehicle 
during Atlas booster phase was essentially zero angle of attack in order to minimize 
aerodynamic heating and loads. Once the vehicle has exited from the sensible atmo- 
sphere, explicit guidance equations were used to steer the vehicle into parking orbit and 
later onto the desired transfer trajectory. The guidance equations used are presented 
in reference 6. 

The direct ascent lunar, Mars, and Jupiter cases used a three-stage version of the 
Atlas-Centaur launch vehicle. For the Jupiter case, the third stage was a Burner n, 
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which is a small solid stage. For the lunar and Mars cases, the third stage was an as- 
sumed cryogenic upper stage. 

In addition to the reference trajectory profiles, a guidance hardware error model 
and error values are required. The error model and error values used in all cases are 
presented in table I. These errors are used to calculate the errors in sensed accelera- 
tion, which in turn, are used as driving functions in the integration of the linearized 
equations of motion. 

The results obtained are summarized in table m. Seven check cases are presented, 
five planetary and two lunar. For the planetary cases, the best agreement is obtained 
for the Mercury case, where the error is 0. 2 percent. The errors for the Venus and 
Jupiter cases are 4 and 6 percent, respectively; the two Mars cases are in error by 
3. 5 and 14 percent. The analytical and N Body results agree within 2. 5 percent for the 
lunar cases. 

The Mars N Body results were obtained for midcourse correction times of 2, 3, 5, 
and 10 days after injection, in order to demonstrate that the required correction velocity 
does not change appreciably for small correction times compared to the total trip time. 
These results are shown in figure 3. 


TABLE m. - SUMMARY OF ERROR ANALYSIS RESULTS 


Case 

Mission 

Midcourse 

correction 

time 

Figure of merit, m/sec 



Miss plus time 

Miss only 




Analytical 

N Body 

Analytical 

N Body 

1 

Mars 

2. 7 Days 

3.90 

3.84 

— 



2 

Mars 

2. 7 Days 

3.39 

2.97 

— 

— 

3 

Jupiter 

1. 1 Days 

8.07 

7.59 

— 

— 

4 

Venus 

2. 6 Days 

3.78 

3.63 

— 

— 

5 

Mercury 

1. 6 Days 

4.81 

4. 80 

— 

— 

6 

Lunar 

20 Hr 

6. 11 

6. 15 

3.57 

3.66 

7 

Lunar 

20 Hr 

6.03 

5.90 

3.66 

3.57 



Figure 3. - Variation of required midcourse 
correction velocity with correction time 
for direct ascent Earth -Mars trajectories. 
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TABLE IV. - COMPARISON OF ANALYTICAL AND N BODY RESULTS 
FOR VENUS TRAJECTORY - CASE 4 
[Figure of merit, 3. 63 m/sec for miss plus time of flight]] 


(a) Analytical results 


Error 

Figure of merit 

Figure of merit 

Percent of 

Numbe] 

r Magnitude 

sensitivity 

coefficient 

per error source 
(miss plus time), 

figure of merit 
squared 



(miss plus time 

1 m/sec 

(miss plus time) 

1 

0.0540 

2.83 

0. 15 

0. 16 

2 

.0540 

12.06 

.65 

2.97 

3 

.0900 

3.37 

.30 

.64 

4 

.0500 

6.50 

.32 

.74 

5 

.0500 

11.44 

.57 

.29 

6 

.0550 

2.64 

. 15 

. 15 

7 

.0420 

14.32 

.60 

2. 53 

8 

.0420 

4.21 

. 18 

.22 

9 

.0420 

11.05 

.46 

1.51 

10 

.0051 

241.34 

1.23 

10.61 

11 

.0051 

15.15 

.08 

.04 

12 

.0051 

102. 10 

.52 

1.90 

13 

.0840 

5.67 

.48 

1.59 

14 

.0930 

20.02 

1.86 

24.28 

15 

.0940 

7.02 

.66 

3.05 

16 

. 1730 

4.93 

.85 

5.09 

17 

. 1900 

4.76 

.90 

5.72 

18 

. 1770 

4.14 

.73 

3.75 

19 

. 1060 

1.18 

. 12 

. 11 

20 

. 1010 

19.72 

1.99 

27.76 

21 

. 1140 

1.68 

. 19 

.26 

22 

.0090 

2.27 

.02 

.00 

23 

.0090 

9.62 

.09 

.05 

24 

.0090 

1.68 

.02 

.00 

28 

.0094 

37. 19 

. 35 

.86 

29 

.0094 

.62 

.01 

.00 

30 

.0094 

21.81 

.20 

.29 

34 

.0130 

6.09 

.08 

.04 

35 

.0130 

.53 

.01 

.00 

36 

.0130 

2.67 

.03 

.01 

37 

.0120 

9. 15 

. 11 

.08 

38 

.0120 

2.48 

.03 

.01 

39 

.0120 

11.27 

. 14 

. 13 

40 

.0240 

22.68 

.54 

2.07 

41 

.0260 

6. 10 

.16 

. 18 

42 

.0290 

11.05 

.32 

.72 

43 

.0090 

2.29 

.02 

.00 

44 

. 0090 

10.07 

.09 

.06 

45 

.0090 

14. 39 

. 13 

. 12 

46 

.0080 

1.77 

.01 

.00 

47 

.0080 

1.82 

. 01 

.00 

48 

.0080 

3.44 

.03 

.0! 
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TABLE IV. - Concluded. COMPARISON OF ANALYTICAL AND N BODY 


RESULTS FOR VENUS TRAJECTORY - CASE 4 
[Figure of merit, 3. 78 m/sec for miss plus time of flightT] 
(b) N Body results 


Error 

Figure of merit 
sensitivity 
coefficient 
(miss plus time) 

Figure of merit 
per error source 
(miss plus time), 
m/sec 

Percent of 
figure of merit 
squared 

(miss plus time) 

Number 

Magnitude 

1 

0.0540 

3. 14 

0. 17 

0.22 

2 

.0540 

11.67 

.63 

3.02 

3 

.0900 

1. 56 

. 14 

. 15 

4 

.0500 

6.00 

.30 

.68 

5 

.0500 

10.93 

. 55 

2.27 

6 

.0550 

2.52 

. 14 

. 15 

7 

.0420 

15. 56 

.65 

3.25 

8 

.0420 

3.77 

. 16 

. 19 

9 

.0420 

10.58 

.44 

1.50 

10 

.0051 

251.84 

1.28 

12.54 

11 

.0051 

13.99 

.07 

.04 

12 

.0051 

105.54 

.54 

2.20 

13 

.0840 

6.61 

.56 

2.34 

14 

.0930 

19.43 

1.81 

24.81 

15 

.0940 

.33 

.03 

.01 

16 

. 1730 

5.52 

.95 

6.93 

17 

. 1900 

4.55 

.86 

5.68 

18 

. 1770 

2.09 

.37 

1.04 

19 

. 1060 

1.32 

. 14 

. 15 

20 

. 1010 

18.87 

1.91 

27.61 

21 

. 1140 

.33 

.04 

.01 

22 

.0090 

2.53 

.02 

.00 

23 

. 0090 

9. 18 

.08 

. 05 

24 

.0090 

.73 

.01 

. 00 

28 

.0094 

38.91 

.37 

1.02 

29 

.0094 

.58 

.01 

.00 

30 

. 0094 

21.24 

.20 

.30 

34 

.0130 

6.56 

.09 

.06 

35 

.0130 

.73 

.01 

.00 

36 

.0130 

2.82 

.04 

.01 

37 

.0120 

9.58 

. 12 

. 10 

38 

.0120 

2.33 

.03 

.01 

39 

.0120 

11.81 

. 14 

. 15 

40 

. 0240 

23.65 

.57 

2. 45 

41 

.0260 

5.62 

. 15 

. 16 

42 

.0290 

10.58 

.31 

.72 

43 

.0090 

2.40 

.02 

.00 

44 

.0090 

9.41 

.08 

.05 

45 

. 0090 

13.65 

. 12 

. 11 

46 

.0080 

1.90 

.02 

.00 

47 

.0080 

2.75 

.02 

.00 

48 

.0080 

3.27 

.03 

.01 
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TABLE V. - COMPARISON OF ANALYTICAL AND N BODY RESULTS FOR LUNAR TRAJECTORY - CASE 6 
[Figure of merit, 3. 57 m/sec for miss only and 6. 11 m/sec for miss plus time of flight^] 


(a) Analytical results 


Error 

Figure of merit 

Figure of merit 

Percent of 
figure of merit 

Figure of merit 

Figure of merit 

Percent of 

Numbe] 

Magnitude 

sensitivity 

per error source 

sensitivity 

per error source 

figure of merit 

coefficient 
(miss only) 

(miss only), 
m/sec 

squared 
(miss only) 

coefficient 
(miss plus time] 

(miss plus time), 
m/sec 

squared 

(miss plus time] 

1 

0.0540 

3.88 

0.21 

0.34 

6.24 

0.34 

0.30 

2 

.0540 

15.38 

.83 

5.41 

28.58 

1. 54 

6.38 

3 

.0900 

1.05 

.09 

.07 

1.11 

. 10 

.03 

4 

.0500 

7. 87 

.39 

1.21 

14.92 

.75 

1.49 

5 

.0500 

12.97 

.65 

3. 30 

22.30 

1. 12 

3.33 

6 

.0550 

2. 81 

. 15 

. 19 

4.84 

.27 

. 19 

7 

.0420 

13. 11 

.55 

2.38 

30. 37 

1.28 

4.36 

8 

.0420 

3. 11 

. 13 

. 13 

6.70 

.28 

.21 

9 

.0420 

10.43 

.44 

1.51 

17.68 

.74 

1.48 

10 

.0051 

318.79 

1.63 

20.73 

680.44 

3.47 

32.28 

11 

.0051 

17.03 

.09 

.06 

32.30 

. 16 

.07 

12 

.0051 

43.59 

.22 

. 39 

95.50 

.49 

.64 

13 

.0840 

3.94 

.33 

.86 

5.36 

.45 

.54 

14 

.0930 

17. 17 

1.60 

19.99 

24.76 

2.30 

14.21 

15 

.0940 

1.26 

. 12 

. 11 

1.29 

.12 

.04 

16 

. 1730 

4.22 

.73 

4. 19 

5.65 

.98 

2.56 

17 

. 1900 

4.02 

.76 

4.58 

5.71 

1.09 

3. 16 

18 

. 1770 

1.33 

.24 

.44 

1.38 

.24 

. 16 

19 

. 1060 

.92 

. 10 

.07 

1.24 

. 13 

.05 

20 

. 1010 

18.47 

1.87 

27.29 

26. 17 

2.64 

18.72 

21 

. 1140 

.30 

.03 

.01 

. 31 

.04 

.00 

22 

.0090 

1.93 

.02 

.00 

2.77 

.02 

.00 

23 

.0090 

8.30 

.07 

.04 

12.74 

. 11 

.04 

24 

.0090 

.42 

.00 

.00 

.43 

.00 

.00 

28 

.0094 

50.64 

.48 

1. 78 

108.47 

1.02 

2. 79 

29 

.0094 

.62 

.01 

.00 

1. 15 

.01 

.00 

30 

.0094 

7.87 

.07 

.04 

12.38 

. 12 

.04 

34 

.0130 

3.94 

.05 

.02 

9.40 

. 12 

.04 

35 

.0130 

. 41 

.01 

.00 

.51 

.01 

.00 

36 

.0130 

1.45 

.02 

.00 

3. 16 

.04 

.00 

37 

.0120 

11. 11 

. 13 

. 14 

23.80 

.29 

.22 

38 

.0120 

2.80 

.03 

.01 

5.22 

.06 

.01 

39 

.0120 

6.56 

.08 

.05 

14.30 

. 17 

.08 

40 

.0240 

27.26 

.65 

3.36 

58.26 

1.40 

5.24 

41 

.0260 

6.80 

. 18 

.24 

12.78 

.33 

.30 

42 

.0290 

10.43 

.30 

.72 

17.68 

.51 

.70 

43 

.0090 

2.45 

.02 

.00 

5.24 

.05 

.01 

44 

.0090 

12.77 

. 11 

. 10 

23.81 

.21 

.12 

45 

.0090 

18.52 

. 17 

.22 

30.57 

.28 

.20 

46 

.0080 

.87 

.01 

.00 

2.07 

.02 

.00 

47 

.0080 

1.87 

.01 

.00 

2.30 

.02 

.00 

48 

. 0080 

4.05 

.03 

.01 

6.67 

.05 

.01 
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TABLE V. - Concluded. COMPARISON OF ANALYTICAL AND N BODY RESULTS FOR LUNAR TRAJECTORY - CASE 6 


[Figure of merit, 3. 66 m/sec for miss only and 6.15 m/sec for miss plus time of flight] 


(b) N Body results 


Error 

Figure of merit 

Figure of merit 

Percent of 

Figure of merit 

Figure of merit 

Percent of 

Number 

Magnitude 

sensitivity 

coefficient 

per error source 
(miss only), 

figure of merit 
squared 

sensitivity 

coefficient 

per error source 
(miss plus time), 

figure of merit 
squared 



(miss only) 

m/sec 

(miss only) 

(miss plus time) 

m/sec 

(miss plus time) 

1 

0.0540 

3.97 

0.21 

0.34 

6.24 

0.34 

0.30 

2 

.0540 

15.91 

.86 

5.51 

28.80 

1.56 

6.39 

3 

.0900 

1.02 

.09 

.06 

1. 11 

. 10 

.03 

4 

.0500 

8.21 

.41 

1.26 

15.19 

.76 

1.53 

5 

.0500 

13.34 

.67 

3.32 

. 22.38 

1. 12 

3.31 

6 

.0550 

2.89 

. 16 

. 19 

4.85 

.27 

. 19 

7 

.0420 

13.94 

.59 

2.56 

31.46 

1.32 

4.61 

8 

.0420 

3.28 

. 14 

. 14 

6.92 

.29 

.22 

9 

.0420 

10.71 

.45 

1.51 

17.70 

.74 

1.46 

10 

.0051 

336.04 

1.71 

21.90 

696.05 

3.55 

33.30 

11 

.0051 

17.76 

.09 

.06 

32.90 

. 17 

.07 

12 

.0051 

46.46 

.24 

.42 

98.77 

.50 

.67 

13 

.0840 

3.96 

.33 

.83 

5.27 

.44 

.52 

14 

.0930 

17.35 

1.61 

19.41 

24.46 

2.27 

13.67 

15 

.0940 

1.23 

. 12 

. 10 

1.28 

. 12 

.04 

16 

. 1730 

4.24 

.73 

4.01 

5.55 

.96 

2.44 

17 

. 1900 

4.06 

.77 

4.43 

5.64 

1.07 

3.03 

18 

. 1770 

1.30 

.23 

.40 

1.37 

.24 

. 16 

19 

. 1060 

.93 

. 10 

.07 

1.21 

. 13 

.04 

20 

.1010 

18.64 

1.88 

26.42 

25.81 

2.61 

17.96 

21 

. 1140 

.30 

.03 

.01 

.31 

.03 

.00 

22 

.0090 

1. 95 

.02 

.00 

2.73 

.02 

.00 

23 

.0090 

8.43 

.08 

.04 

12.65 

. 11 

.03 

24 

.0090 

.41 

.00 

.00 

. 43 

.00 

.00 

28 

.0094 

53.38 

. 50 

1.88 

110.96 

1.04 

2. 87 

29 

.0094 

.64 

.01 

.00 

1. 17 

.01 

.00 

30 

.0094 

8.01 

.08 

.04 

12.30 

. 12 

.04 

34 

.0130 

4. 15 

.05 

.02 

9.60 

. 12 

.04 

35 

.0130 

.41 

.01 

.00 

. 50 

.01 

.00 

36 

.0130 

1. 54 

.02 

.00 | 

3.26 

.04 

.00 

37 

.0120 

11.71 

. 14 

. 15 1 

24. 34 

.29 

.23 

38 

.0120 

2.92 

.04 

.01 

5.31 

.06 

.01 

39 

.0120 

6.99 

.08 

.05 

14.78 

. 18 

.08 

40 i 

.0240 

28.73 

.69 

3.55 

59.60 

1.43 

5.41 

41 

.0260 

7.08 

. 18 

.25 

13.01 

.34 

.30 

42 

.0290 

10.71 

.31 

.72 

17.70 

.51 

.70 

43 

.0090 

2. 58 

.02 

.00 

5.36 

.05 

.01 

44 

.0090 

13.29 

. 12 

. 11 

24.22 

.22 

. 13 

45 

.0090 

18.95 

. 17 

.22 

30. 53 

.27 

.20 

46 

.0080 

.91 

.01 

.00 

2.11 

.02 

.00 

47 

.0080 

1.87 

.01 

.00 

2.28 

.02 

.00 

48 

.0080 

4. 14 

.03 

.01 

6.66 

.05 

.01 
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TABLE VI. - 

FIGURE OF MERIT COVARIANCE MATRIX (MISS PLUS TIME) FOR VENUS TRAJECTORY 

(a) Analytical partials 

- CASE 4 

0. 7345719x10’^ 
.5974771X10' 6 
1421085X10' 12 
. 1587145X10’ 2 
. 8998715X10" 2 
7227641X10" 10 

0. 5974771X10’ 6 
. 199742 5xl0“ 6 
-.2664535X10' 14 
. 3208562X10" 3 
.6506040X10' 3 
1032960X10" 9 

-0. 1421085X10’ 12 
2664535X10’ 14 
. 1997425X10' 6 
1046427X10" 7 
1594759X10' 6 
.6915269X10' 4 

0. 1587145X10' 2 
.3208562X10' 3 
1046427X10' 7 
.5862252 
1.8413685 
-.1117587X10' 7 

0.8998715X10" 2 
. 6506040X10' 3 
1594759X10' 6 
1. 8413685 
11.067776 
-.2384186X10' 6 

-0.7227641X10* 10 
-.1032960X10' 9 
.6915269X10' 4 
-.1117587X10' 7 
-. 2384186X10” 6 
. 2394136X10* 1 

(b) N Body partials 

0. 7456713x10’^ 

0.665U50X10' 6 

-0.6675128X10' 7 

0.1764756X10' 2 

0.9074949X10' 2 

0. 3023507X10' 4 

. 6651150x10”® 

.2207393X10' 6 

1500603X10' 7 

.3617413X10' 3 

. 7209460X10” 3 

. 8065659X10' 5 

-.6675128X10' 7 

1500603X10" 7 

. 2069464X10' 6 

-.2712879X10' 4 

-. 7663904X10” 4 

-. 1109500X10" 3 

. 1764756X10' 2 

.3617413X10' 3 

-.2712879X10' 4 

.6763185 

2.0356938 

.14997 15X10' 1 

. 9074949X10' 2 

.7209460X10" 3 

-.7663904X10' 4 

2.0356938 

11.092905 

. 4544180X10' 1 

.3923507X10' 4 

.8065659X10' 5 

1109500X10" 3 

. 1499715X10' 1 

.4544180X10' 1 

.5948564X10' 1 





TABLE VII. - FIGURE OF MERIT COVARIANCE MATRIX FOR LUNAR TRAJECTORY - CASE 6 


(a) Analytical partials; miss plus time 

0. 4674652X10’ 4 

0. 1270308X10' 5 

-0. 6821210X10' 12 

0. 1515782X10' 2 

0. 5531167X10' 1 

-0. 1564331X10' 9 

. 1270308X10' 5 

. 1648680X10' 6 

-.7105427X10' 14 

. 1958213X10' 3 

. 1488603X10' 2 

-.4684253x10' 11 

6821210X10’ 12 

-.7105427x10' 14 

. 1543766X10' 6 

-.3613799X10" 8 

5540254X10' 6 

-.1602241X10' 4 

. 1515782X10" 2 

. 1958213X10' 3 

-.3613799X10" 8 

.2325870 

1.7761296 

-.4656613X10' 8 

. 5531167X10’ 1 

. 1488603X10' 2 

5540254X10' 6 

1.7761296 

65. 448912 

-.9536743X10' 6 

1564331X10' 9 

-.4684253X10" 11 

1602241X10' 4 

-.4656613X10" 8 

-. 9536743X10' 6 

. 1662929X10' 2 

(b) Analytical partials; miss only 

0. 1131277X10' 4 

-0.2109825X10' 6 

0. 1002750X10" 5 

-0.2481504X10' 3 

0. 134440 1X10' 1 

-0. 1040733X10' 3 

. 2109825X10’ 6 

. 1029433X10' 6 

. 4191951x10"^ 

. 1220812X10' 3 

-.2620735X10' 3 

-.4350733X10" 5 

. 1002750X10' 5 

. 4191951x10’^ 

. 1259994X10' 6 

.4991549X10' 4 

. 1184639X10' 2 

-. 1307721X10" 4 

-.2481504X10" 3 

. 1220812X10' 3 

. 4991549X10' 4 

. 1447765 

-.3081134 

-. 5180882X10’ 2 

. 1344401X10' 1 

-.2620735X10' 3 

. U84639X10' 2 

-.3081134 

15.977962 

-. 1229724 

1040733X10" 3 

4350733X10' 5 

-.1307721X10" 4 

5180882X10' 2 

-. 1229724 

. 1357254X10" 2 



(c) N Body partials; 

miss plus time 



0. 4869086X10’ 4 

0. 1181656X10' 5 

-0. 1495929X10' 6 

0. 1791531X10' 2 

0. 5762047X10' 1 

0. 1037888X10" 4 

. 1181656xl0' 5 

. 1584029X10' 6 

5109898X10' 8 

. 1973609X10' 3 

. 1382835X10' 2 

. 2036971X10" 6 

-.1495929X10" 6 

5109898X10' 8 

. 1550908X10' 6 

-.7269783X10" 5 

-. 1770799X10' 3 

-. 1587410X10’ 4 

.1791531X10" 2 

. 1973609X10' 3 

-.7269783X10" 5 

.2484564 

2. 1016650 

. 3258169X10" 3 

. 5762047X10" 1 

. 1382835X10' 2 

1770799X10" 3 

2. 1016650 

68. 189584 

.1231161X10" 1 

. 1037888X10" 4 

. 203697 lxlO' 6 

1587410X10' 4 

.3258169X10" 3 

. 1231161X10' 1 

. 1625580X10' 2 


(d) N Body partials; miss only 


0. 1247782X10' 4 

-0.2527232X10' 6 

0. 1000962X10" 5 

-0. 1999562X10' 3 

0. 1483115X10' 1 

-0. 1029777X10' 3 

-.2527232X10' 6 

. 101587 9xlO' 6 

. 4046299X10' 7 

. 1184792X10' 3 

-. 3120274X10" 3 

-.4286299X10' 5 

. 1000962X10' 5 

. 4046299X10' 7 

. 1185355X10’ 6 

. 5600343X10' 4 

. 1182416X10" 2 

-. 1227255X10" 4 

-. 1999562X10' 3 

. 1184792X10' 3 

. 5600343X10" 4 

. 1389373 

-.2514756 

-.5908076X10' 2 

.1483115X10' 1 

-. 3120274X10" 3 

. 1182416X10' 2 

-.2514756 

17.629724 

-.1216305 

-. 1029777X10' 3 

-.4286299X10' 5 

-. 1227255X10" 4 

-. 5908076X10' 2 

-.1216305 

. 1270743X10' 2 
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to 
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The FOM results presented in table in are obtained by statistically combining the 
effects of the guidance error sources presented in table I. The contribution of each error 
source to the total FOM is presented for Venus and lunar missions (cases 4 and 6) in 
tables IV and V. In addition to the FOM per error source, the FOM sensitivity coeffi- 
cient (FOM per unit error source) and percent of FOM squared are also presented for 
each error source. These quantities allow the FOM to be quickly recalculated for differ- 
ent error values. The analytical and N Body FOM covariance matrices are also pre- 
sented for the same two cases in tables VI and VII. 


Sample Problem 

Suppose it is desired to determine the trajectory resulting in the largest FOM for a 
particular planetary launch opportunity. Such a problem arises in preliminary design 
studies, where it is necessary to size the midcourse propulsion system. If the analytical 
equations described in this report are used, the possible trajectories for the opportunity 
may be described in terms of the powered flight trajectory variables. For a two-burn 
mission, these variables are launch azimuth, energy, and parking orbit coast time. For 
a direct ascent mission, injection true anomaly replaces coast time as a variable. The 
largest FOM can then be determined by searching all possible values of the powered 
flight variables, the range of which may be found in the literature; for example, refer- 
ences 4 and 5. 

The use of this technique is illustrated on figure 4 for a two-burn Mars mission. 
Figure of merit is plotted against injection energy for different values of parking orbit 
coast time. The launch azimuth is constant on figure 4; therefore, other figures must 
be generated for different launch azimuths. These figures are not presented here. The 
largest FOM may be easily determined from the figures. 



Figure 4. - Variation of figure of merit with energy and coast time. 
Launch azimuth, constant; two-burn Mars mission. 
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After the values of the powered flight variables resulting in the largest FOM have 
been determined, an N Body trajectory corresponding to this case can be generated, if 
a more accurate estimate of FOM is desired. Even if the analytical FOM is in error by 
a few percent, the trends resulting in the choice of the worst case powered flight vari- 
ables can be assumed to be correct. 


CONCLUSIONS 

Analytical equations have been developed for computing partial derivatives relating 
midcourse correction velocity requirements to injection errors for planetary and lunar 
missions. These partial derivatives are used to calculate the guidance system FOM. 

The equations are easily programmed on a digital computer and require very little com- 
puter time. In addition, no free flight reference trajectories are required. Thus, error 
analysis results can be obtained without first having to calculate detailed targeted trajec- 
tories, a process which requires a powered flight -free flight (N Body) targeting program. 

Analytical and N Body results have been generated for two lunar and five planetary 
missions. These results show good agreement between the analytical and N Body results. 
For the planetary cases, the best agreement is obtained for the Mercury case where the 
error is 0.2 percent. The errors for the other planetary cases vary between 2 and 
14 percent. For the lunar cases, the results agree within 2. 5 percent. 

Because of the accuracy of the results presented and the ease with which these re- 
sults may be obtained, the analytical equations provide a useful preliminary design tool 
for estimating space mission FOM. The analytical technique can also be readily ex- 
tended to apply to other mission types. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 12, 1968, 

125-17-05-01-22. 



APPENDIX A 


A,B,] 

c,Dj 

> Danby's matrices 

a 

thrust acceleration vector, 
m/sec^ 

c 

cos 7) (appendix B) 

E 

2 2 

energy, m /sec 

E[ J 

expected value operator 

e 

eccentricity 

e k 

k th guidance hardware error 
source 

G ik 

matrix of partial derivatives, 
Sv i /9S k 

h 

2 

angular momentum, m /sec 

l,m; 

P,Q, 

j, Danby's matrices 

P ik 

matrix of partial derivatives, 
SS i /0e k 

P 

semilatus rectum, m 

r 

radius, m 

r, e,h 

unit radial, tangential, normal 
coordinate system 

S 

state variable 

s 

sin t? (appendix B) 

t 

time 

u, V 

vectors defined in equations (8) 
and (10) 

V 

velocity, m/sec 


magnitude of midcourse correc 
tion velocity, m/sec 

midcourse correction velocity, 
m/sec 

z unit vector pointing at North 

Pole 

a, /3, p coefficients used in appendix C 

y flight path angle, rad 

A( ) error quantity 

5( ) linearized error quantity 

7 ] true anomaly, rad 

A FOM covariance matrix 

u central body gravitational con- 

3 2 

stant, m /sec 

v variable defined in appendix B, 

(sec)"* 

a injection covariance matrix 

cp latitude of moon at arrival with 

respect to earth's equator, 
rad 

\p inclination of lunar orbit to 

earth's equator, rad 

Subscripts: 

A arrival 

a asymptote 

act actual 

c midcourse correction 

e earth escape 


SYMBOLS 

Av, 


Av, 
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h 

normal 

r 

radial 

I 

injection 

t 

tangential 

M 

miss only 

Superscripts: 

MT 

miss plus time 

T 

matrix transpose 

m 

moon 

— 

vector 

mes 

measured by the guidance system 


unit vector 

nom 

nominal 

General: 

P 

perigee 

• 

vector dot product 


x vector cross product 



APPENDIX B 


DANBY'S METHOD AND MATRICES 

Errors in position and velocity at any time t on a coast ellipse are related to errors 
at periapse by the following matrix: 


'5r\ / L(t, t p ) M(t,t p )\ / 5r 


v 5v/ t \P(t,t ) Q(t. tp)/ \5v/ t 


(Bl) 


where L, M, P, and Q are 3x3 matrices, the elements of which are presented later 
in this appendix as functions of e, r p , r, r), p, and (t - t p ). To relate times t^ and t 2 
(neither at periapse), the following relation is applied: 


'6r 


/A(t 2 , tj) B(t 2 , tjA M 


v6v/ t \C(t 2 , tj^) Dft^tj)/ \8v /t 
2 1 


(B2) 


where A, B, C, and D are 3x3 matrices and are of the form: 


-\ 


A^’tf) - LCtgjtp) 

Q T (t 1 ,t p ) - M(t 2 ,t p ) 

P (tptp) 

Bft^tp = M(t 2 ,t p ) 

L T <t r , t p ) - L(t 2 ,t p ) 

M T (t 1 ,t p ) 

C^tj) = p(t 2 , t p ) 

Q T (tj,t p ) -Q(t 2 ,t p ) 

P T (ti,t p ) 

Dftg.tj) = Q(t 2 ,t p ) 

L T (‘l>‘p> -P(t 2 ,t p ) 

M T (t 1 ,t p ) 


(B3) 


The elements of the L, M, P, and Q matrices are given by 

Zv 


J 11 


r p (l - e) 


'c 2 + c(2 - e) - 2] + 


— — S(t - t ) 
(1 - e) (1 + e) 1 / 2 
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1 1 III II III III I II 



L 12 = 


r p (l + e) 


s(l - c) 


J 21 


V 1 - e > 


s(c + 2) - 


Zv 


— (c + e) (t - t ) 
(1 - e) (1 + e) 1 / 2 


L 22 “ 


r p(l + e) 


[c 2 - c(l - e) + 1 


L 33 c 

r p 


M 


11 


^r p (l + e) 3/2 


s(-c + 2 + e) 


M 


12 


w p (l - e) (e + 1) 


2(c - 1) (c + 2) + s(t - t ) 

1/2 1 - e p 


M 


21 


r (1 - c) 2 


^r p (l + e) 3 / 2 


M 22 " 


w (1 - e) (1 + e) 


2s(c + 1 + e) (c + e) (t - t ) 

1/2 1 - e p 


M 33 " 


^r p (e + l) 1 / 2 


3^ 2 r 2 


P ll = 


(1 - e) (e + l) 1 / 2 


s(-ec^ - 2c + 1 - e) + 2— c(t - t ) 

r 2 (l - e) P 


12 


(1 + e ) 3 / 2 


(-ec 2 - 2c 2 + c + 1 + e) 
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21 


(1 - e) (e + 1) 


1/2 


3i^r^ 

(ec 3 + 2c 2 - c - 1 - e) + 2_ (t - t ) 

r 2 (l - e) P 


P 00 = t s(-ec 2 - 2c + 1) 


22 


(1 + e) 3/2 


P 33 " 


-v 


(1 + e) 


1/2 


Qjj = — [^-ec 3 - 2c 2 + 2c + ec + (1 + e) 2 J 

(1 + e) 2 


2>vy £ 


v 1/2 


Qio = — s(-2ec 2 - 4c + 1 - e) + — "~ p (1 + e) A/ * c(t - t ) 

v2 2 /- v P 


(1 - e) 


r (1 - e) 


Q 21 = 7 s(l - c) (ec + 2 + e) 


(1 + e)^ 


3w' 


Q 22 = — j"2ec 3 + 4c 2 - c(l + e) -2 - e - e 2 J + 2 — (l + e J^ 2 s(t - t ) 

(1 - e 2 ) r 2 (l - e) 


Qoq (c + e) 

33 1 + e 

where 

s = sin 7] 

C = COS 7j 

v = 

All matrix elements not listed are equal to zero. 
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APPENDIX C 


ANALYTICAL EQUATIONS FOR CALCULATION OF FOM COVARIANCE MATRIX 

Planetary case: 

Calculate reference trajectory variables 


h = r x v 



M 


E = - — 

2 r 


V= 


e = */l + 


2Ep 


r 7 a = cos 


•1 (- 
\ e/ 




r » v 
er 


r = — 


1 + e 


l I " *P ‘ 


-i/Z / 1 \ / i 

V M \l - e 2 / Ve 2 - 1 


^Ve + 1 + V e - 1 tan — 


log 


^V e + 1 - Ve - 1 tan — 

2 , 


t - t = 0 
a p 


e sin r/j 


1 + e cos ? 7 j 


(Cl) 


Calculate Danby’s matrices L Jp , M Jp , Pj p , Q Ip , L ap , M &p , P ap , Q ap using the last 
equation in appendix B. 
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Calculate radial coordinate system at perigee. 


v x h _ r_ 
ep. er 


r x v 
Ir x v| 


W r p 


Calculate Danby's matrices C a j and D a j using equation (B3) in appendix B. 


G = < c al D al> 


A = G T G 


The injection error vector must be expressed in the perigee coordinate system 


5r • r 


5r • h. 


6v • r 


5v • 8 


6v • h 


Lunar case: 


Calculate reference trajectory variables using equations (Cl) except that the equa- 
tions for 77 , (t T - t ), and (t - t ) are replaced by 

dip dp 


-1 "m 


1,1 ■ ill mu i ii mi i 
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*A " same as equation (C 5 ) but with 77 A instead of ?jj 

t c - tp same as equation (C 5 ) but with r) c instead of 7]j 

Iterate on rj c to satisfy t c - tj equals the desired midcourse correction time. 

Calculate Danby's matrices L, p , M Ip , P, p , Q Ip , L cp , M cp , P cp , Q^, L mp , M mp> 

^mp’ am * ^rap' 

Calculate radial coordinate system at perigee using equations (C 2 ). Calculate posi- 
tion and velocity at lunar arrival. 


r A = 


1 + e cos 




(cos ?7 A r p + sin V A $ p ) 


e - h x r A 

v = _ h x r + 

P P pr A 


Calculate lunar velocity vector 


'■■ft 


-cos 
cos <p 


r x z ± 
m 




C0S * ?m * <?m >< 


COS ^(p 


where 


sin cp = r • z 
^ m 


cos (p = Vl - sin^tp 

Express (v A - v ) in radial-perigee coordinate system. 


&A ' ? m> ' f p 
^A - ? m> = ^A - • d p 

&A - ' % 
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Calculate Danby's matrices A mI , B mI , and B mr , using equations (B3) in appendix B 


me 


Miss plus time correction: 


G = -' B mc A mI 




Miss only correction: 


A MT “ G G 


V = -B _1 (v A - v ) 
mc v A nr 


a m " a mt 


V T G 
’ V T V 


The injection error vector is expressed in the perigee -radial coordinate system, as in 
equation (C3). 
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APPENDIX D 


CALCULATION OF UNIT LUNAR VELOCITY VECTOR, v 


m 


Let v m be expressed as 


% = “ r m + ^ r m X z > + P r m X (f m X z > < D1 ) 


Three equations are available for determining a, /3, and p based on the assumption that 
the moon's orbit about the earth is circular with an inclination if/ to the earth's equator. 


v v =1 
m m 


v • r =0 
m m 


(? m X 2 = -*m ‘ (f m X z > = cos ^ 


Combining equations (Dl) and (D2) results in 


(D2) 

(D3) 

(D4) 


+ ^m X z > 2 + P 2 [ f m X (? m X = 1 


~\l2 


(D5) 


since the three vectors used in defining v m are mutually orthogonal. But 

(? m x z) 2 = 1 - (? m • z) 2 = 1 - sin 2 cp = cos 2 <p 
where <p is the latitude of the moon relative to the earth. Also, 

Pm X tfm x2) ] 2 = ( f m X z > 2 = cos2( P 


Therefore, equation (D5) becomes 


a 2 + ()S 2 + p 2 )cos 2 cp = 1 


From equations (Dl) and (D3) it follows that 


r • = (?„ • r )a = 0 

m m m m' 


which implies that 


(D6) 


(D7) 
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a = 0 


(D8) 


From equations (Dl) and (D4) 

v m ' X ^ = X ^ = _cos ^ 


or 


(3 =■ (D9) 

COS 2 <£> 

where equation (D6) has been used. Combining equations (D7) and (D9) results in 


P = 


±1 

cos cp 



COS^l// 

cos 2 #> 


(DIO) 


Substituting equations (D8) to (DIO) into equation (Dl) results in 


v m 


cos cp 




cos cp 


cos 2 <p 


(Dll) 


The ambiguity in sign in equation (Dll) indicates that there exist two possible velocity 
vectors satisfying equations (D2) to (D4). These two possibilities correspond to a rising 
or setting moon. Taking the dot product of equation (Dll) with z results in 


m 


z = ± — — — l/l - [(? • z)r - z] ■ z = =f Vcos 2 

cos <P V cos 2 <p 


cp - COS*V 


(D12) 


Therefore, the sign choices in equation (Dll) correspond to 
+ sign indicates descending moon 

- sign indicates ascending moon 
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